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ABSTRACT 

We report the discovery of two giant planets orbiting stars in Praesepe (also known as the Beehive 
Cluster). These are the first known hot Jupiters in an open cluster and the only planets known to 
orbit Sun-like, main- sequence stars in a cluster. The planets are detected from Doppler shifted radial 
velocities; line bisector spans and activity indices show no correlation with orbital phase, confirming 
the variations are caused by planetary companions. Pr0201b orbits a V — 10.52 late F dwarf with 
a period of 4.4264 ± 0.0070 days and has a minimum mass of 0.540 ± 0.039 Mjup, and Pr0211b 
orbits a. V = 12.06 late G dwarf with a period of 2.1451 ± 0.0012 days and has a minimum mass 
of 1.844 ± 0.064 Mjup. The detection of 2 planets among 53 single members surveyed establishes 
a lower limit on the hot Jupiter frequency of 3.8^2 4 in this metal-rich open cluster. Given the 
precisely known age of the cluster, this discovery also demonstrates that, in at least 2 cases, giant 
planet migration occurred within 600 Myr after formation. As we endeavor to learn more about the 
frequency and formation history of planets, environments with well-determined properties - such as 
open clusters like Praesepe - may provide essential clues to this end. 

Subject headings: open clusters and associations: individual (Praesepe, M44, NGC 2632, Beehive) — 
planetary systems — stars: individual (BD-h20 2184, 2MASS J08421149-M916373) 



1. INTRODUCTION 

Exoplanet studies over the last 15 years have demon- 
strated that at least 10% of FGK stars harbor gas gi- 
ant planets, with many of them at surprisingly small 
separations, im plying inward migration after formation 
([Wright et aLl uOll). Although the mechanism by which 
most planets migrate is not understood, powerful con- 
straints on proposed theories of migration can be estab- 
lished by determining the orbital properties of planets 
at young or adolescent ages (< 1 Gyr). For example, 
if migration occurs prima rily due to interactions with 
a circumstellar disk (e.g., IGoldreich fc Tremaind [19801 : 
ILin et all I1996[) . the migr ation must occur befo re the 
disk dissipates (~10 Mvr: [Carpenter et al.l 12006*). and 
is predicted to circularize orbits. Alternatively, if mi- 
gratio n occurs primarily due to planet-planet scattering 
(e.g., I Adams fc LaughlinI I2003D . the process may take 
hundreds of millions of years to occur and can produce 
highly eccentr ic orbits, prior to an y tidal circularization 
(see review bv lLubow fc"ldal[20T0l) . 

A direct way to find planets that can potentially be 
used to constrain theories of migration is to search for 
them in young open clusters. However, until now, only 2 
open cluster stars were known to harbor planets - e Tau 
in the Hyades (Sato et al. 2007) and TYC 5409-2156-1 
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in NGC 2423 (|Lovis fc Mavoii[2007[) - both of which are 
giant stars and thus, by necessity, have planets on wider 
orbits than those occupied by hot Jupiters. In both 
cases the host stars are of intermediate mass (2.7 and 
2.4 Mq), likely A or B type stars when on the main se- 
quence. The lack of detected planets orbiting FGK main 
sequence stars (which are often referred to as Sun-like 
stars) in open clusters has remained despite radial veloc- 
ity (RV) suvey s of 94 dwarfs in the metal-rich Hyades 
([Paulson et a l.l [200I mean Fe/Hl = -HO. 13) and 58 
dwarfs in M67 ([Pasauini et al.l [20121) . as well as numer- 
ous t r ansit searches in other clusters (e.g.JHartnian et al.l 
[200I [Pepper et al.|[20Ql iMocheiska et al.ll2006D . While 
the lack of detec tions may still be the result of small sam- 
ple sizes fe.g.. [van Saders fc Gaudil I2011D . milhmeter- 
wave studies of disks around stars in the Orion star form- 
ing region, which may evolve into an ope n cluster, of- 
fer a plausible astrophysical explanation. [Eisner et al.l 
([2008D find that most solar-type stars in this region do 
not possess disks massive enough to form gas giant plan- 
ets. One may also speculate that for the few stars 
capable of forming planets, the remaining disk masses 
may be insufficient to p ermit inward migration (see also 
iDebes fc JacksM3[20ia ). 

In an attempt to (1) more confidently determine 
whether planet formation and/or migration is inhibited 
around stars within clusters, and (2) potentially discover 
planets with known ages and measurable orbital prop- 
erties, we have carried out an RV survey of stars in 
the Praesepe open cluster. Here we present the semi- 
nal result of that survey - the discovery of the first 2 hot 
Jupiters orbiting Sun-like stars in a cluster. 

2. SAMPLE SELECTION 

Stars were selected from the Praesepe open cluster 
because it is relatively nearby (170 pc), has '--^lOOO 
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TABLE 1 

Target List and Observations Summary 



Star 


a 


<5 


V 


N 


<^obs 




(J2000) 


(J2000) 


(mag) 




(ms-1) 


Pr0044 


08 : 34 : 59.6 


+21 : 05 ; 49.2 


11.06 


6 


42.0 


Pr0047 


08 : 35 : 17.8 


+19 : 38 : 10.2 


12.24 


5 


13.4 


PrOOSl 


08 : 35 : 54.5 


+18 : 08 : 57.8 


10.88 


5 


18.3 



Note. — Table [T] is presented in its entirety in the online 
journal. A portion is presented here for guidance regarding its 
form and content. 

known members, a well determined age (600 Myr; 
Hamblv ct al. 1995: Kraus fc Hillenbran3l2007l:lAn et alj 
2007: Caspar et al.,,2009.: .Delorme et al.N20lH) . and sig- 



nific antly elevated m etallicity ([Fe/H] = +0. 27 ± 0.10 
dexJPace et al.l[200l [Fe/H] = +0.11 ± 0.03, lAn et al.l 
I2007IV Its high metallicity is important because giant 
planet freq uency is strongly correlated with host star 
metallicity (iSantos et al.l 120041: iFischer fc Valentil 120051: 
iJohnson et al.ll2010D : a metallicity as high as +0.27 dex 
implies an increase in the giant planet frequency of a 
factor of nearly 4 relative to solar metallicity. If this cor- 
relation applies to open cluster stars, as many as 1 in 20 
Praesepe stars could harbor a hot Jupiter, and 1 in 400 
could host a transiting giant planet. 

Cluster members w ere selected from the mem ber- 
ship list assembled by iKraus fc Hillenbrandl (|2007[ ) , ex- 
cluding stars with known spectros copic or visua l com- 
panions {M crmilliod ct al. 200^ IBouvier etHI 120011: 
IPatienceet al.l 120021) . To ensure the velocity preci- 
sion would be sufficient to detect substellar companions, 
we limited our initial search to slowly rotating, bright 
stars (wsini < 12 kms~^; V < 12.3). After apply- 
ing these cuts, our sample contained 65 stars. Initial 
RV measurements revealed 12 stars to be obvious spec- 
troscopic binaries (ARV >> 1 kms"'^) or non- members 
(|RV — RVciusterl > 5 kms"-^), leaving 53 viable targets 
for our search. These 53 stars are listed in Table [TJ 

3. OBSERVATIONS 

We us ed the Tilling hast Reflector Echelle Spectrograph 
(TRES: lF{iris3[2008^ mounted on the 1.5-m Tilhnghast 
Reflector at the Fred L. Whipple Observatory on Mt. 
Hopkins, AZ to obtain high resolution spectra of Prae- 
sepe stars, between UT 6-Jan-2012 and 16-Apr-2012. 
TRES is a temperature-controlled, fiber-fed instrument 
with a resolving power of R'-^44, 000 and a wavelength 
coverage of ^-^3850-9100 A, spanning 51 echelle orders. 

We aimed to observe each target on two to three con- 
secutive nights, followed by another two to three consec- 
utive nights ^1 week later. This strategy should be sensi- 
tive to most massive planets with periods up to 10 days. 
Though we were sometimes forced to deviate from the 
planned observing cadence because of weather and in- 
strument availability, we were able to obtain 5-6 spectra 
of each of our 53 targets. Exposure times ranged from 
3-30 minutes, yielding a typical SNR per resolution ele- 
ment of ^40. We also obtained nightly observations of 
the lAU RV standard star HD 65583, which is --14 de- 
grees from Praesepe. Precise wavelength calibration was 
estaiblished by obtaining ThAr emission-line spectra be- 
fore and after each spectrum, through the same fiber as 
the science exposures. 



4. ANALYSIS 

4.1. Spectroscopic Reduction and Cross Correlation 

Spectra were optimally extracted, rectified to intensity 
vs. wavelength, and for each star the individual spectra 
were cross-correlated, order by order, using the strongest 
exposure of that sta r as a template (for details, see 
IBuchhave et al .112 01 01) . We typically used ~25 orders, re- 
jecting those plagued by telluric absorption, fringing far 
to the red, and low SNR far to the blue. For each epoch, 
the cross correlation functions (CCFs) from all orders 
were added and fit with a Caussian to determine the rel- 
ative RV for that epoch. Internal error estimates (which 
include, but may not be limited to, photon noise) for 
each observation were calculated as Oint = RMS(w) / y/N , 
where v is the RV of each order, TV is the number of or- 
ders, and RMS denotes the root-mean-squared velocity 
difference from the mean. 

To evaluate the significance of any potential velocity 
variation, we compared the observed velocity dispersions 
{(Jobs), illustrated in Figure [1] to the combined measure- 
ment uncertainties, which we assumed stem from three 
sources: (1) internal error, amt (described above), (2) 
night-to-night instrumental error, cttreS: ^-nd (3) RV jit- 
ter induced by stellar activity, cr,. 

Before assessing the instrumental error, we first used 
observations of HD 65583 to correct for any systematic 
velocity shifts between runs (such as a 25 ms~^ offset 
caused by a shutter upgrade in mid-March) . This was 
done by determining the median RV of HD 65583 for 
each run, and adjusting all RVs from that run by the 
amount required to make the median RV of HD 65583 
constant over all runs. After applying this correction, the 
RMS of the HD 65583 RVs was 10.8 ms'S with internal 
errors of only 6 ms^^. Since we expect negligibile stellar 
jitter for the RV standard, we estimate the instrumental 
floor error to be cttres = VlO.82 — & ms~^ = 9 ms~^. 

In many cases the observed velocity dispersions are too 
large to be explained by internal and instrumental errors 
alone, implying substantial stellar jitter. We calculate 
the stellar jitter, cr* = V ^obs ~ '^'int ~ ^tres ■ ^^'^ mean 
stellar jitt er is 13 ms~^, which is similar to that found 
by Paulso n" et al.l (16 ms~^) for coeval Hyades members. 

Accounting for internal, instrumental, and stellar 
noise, we constructed a fit of each star's RVs as- 
suming a constant velocity, and then calculated f (x^)i 
the probability that the observed value would arise 
from a star of constant RV. Pr0201 (BD+20 2184) and 
Pr0211 (2MASS J08421149-I-1916373) stood out, with 
P(x^) < 0.001. We obtained additional spectra of these 
stars, and in both cases a Lomb-Scargle periodogram re- 
vealed significant periodicity. Their radial velocities are 
presented in Table [21 

4.2. Orbital Solutions 

We used a Markov Chain Monte Carlo (MCMC) anal- 
ysis to fit Keplerian orbits to the radial velocity data 
of Pr0201 and Pr0211, fitting for orbital period P, time 
of conjunction Tc, the radial velocity semi-amplitude K, 
the center-of-mass velocity 7, and the orthogonal quanti- 
ties Y^cosa; and y^sino;, where e is eccentricity and uj 
is the argument of periastron. We calculated errors from 
the extent of the central 68.3% interval of the MCMC 
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Fig. 1. — Observed velocity dispersions of the 53 stars in our 
sample. The two planet hosts are indicated by the solid red boxes. 
We note that while Pr0201 resides in the upper tail of the distribu- 
tion, other stars with similar RMS values do not necessarily host 
planets; large internal errors and larger than average jitter are two 
possible reasons for a large RMS. For reference, a 1 Afjup planet 
in a 3 day orbit around a 1 Mq star has a full orbital amplitude 
of 281.7 ms-^ 



TABLE 2 

Relative Radial Velocities of Planet Hosts 



BJD 


RV 


O'RV 


BJD 


RV 




(-2455900) 


(m s" 




(-2455900) 


(m s" 




Pr0201 


32.841253 


-34.7 


18.6 


94.824805 


19.0 


20.2 


33.821938 


67.3 


23.9 


96.733645 


139.4 


18.7 


34.857994 


48.0 


21.8 


97.744411 


-66.6 


17.7 


39.945964 


-42.6 


20.2 


98.649923 


-63.6 


29.6 


41.020182 


-49.5 


27.0 


121.701500 


-7.3 


22.9 


57.905296 


-8.2 


21.4 


122.646374 


81.9 


17.3 


58.785164 


-29.3 


20.1 


123.637192 


49.5 


16.1 


59.818953 


0.0 


12.0 


124.646740 


-6.3 


14.0 


60.759119 


30.1 


17.0 


125.693014 


-4.1 


16.4 


70.971372 


-18.7 


24.6 


126.718364 


47.8 


30.4 


71.796286 


-51.7 


42.5 


128.683648 


7.1 


18.3 


72.970059 


13.5 


29.8 


129.690867 


13.7 


18.4 


81.771272 


8.6 


23.5 


130.724779 


89.9 


23.6 


82.697156 


87.5 


15.9 


131.705049 


85.2 


20.7 


84.894042 


-90.5 


18.9 


132.735355 


48.4 


24.8 


86.805010 


27.0 


19.6 


133.638578 


-18.2 


12.0 


87.644768 


102.1 


19.6 








Pr0211 


86.772462 


-75.7 


19.1 


124.663763 


212.0 


21.4 


87.720167 


413.8 


13.8 


125.714331 


113.1 


16.9 


88.721641 


-193.0 


21.3 


126.737245 


275.7 


30.0 


96.712491 


311.8 


24.9 


128.698449 


428.7 


18.2 


97.756832 


0.0 


13.8 


129.716113 


-52.5 


21.0 


98.662801 


405.8 


14.7 


130.741942 


480.6 


26.5 


121.718407 


346.4 


18.9 


131.720488 


-156.9 


50.5 


122.663113 


108.0 


20.3 


132.715068 


504.4 


24.7 


123.652317 


208.1 


21.8 


133.654723 


-98.1 


19.9 


Note. — The errors 


listed here are internal error estimates, 



but in the orbital solutions we include an assumed stellar jitter 
of 13 ms~^ and an instrumental floor error of 9 ms~^, added 
in quadrature with the internal errors. 

posterior distributions. 

The full orbital solutions give eccentricities of e = 
0.156lJ^;?i2 for Pr0201 and e = 0.046t^;^24 for Pr0211. 
However, it can take many precise observations to ac- 
curately measu r e sm all, non-zero eccentricities (e.g., 
iZakamska et al.l 120111 ) . and both are consistent with 
e = to within 2-(T, so we advise caution to not over- 
interpret these results; for short period planets such as 
these, we expect that in the absence of additional bodies, 
tidal forces should have already circularized the orbits 



(e.g.. lAdams fc LaughlinI [2006D . We also note that the 
other orbital parameters are changed by less than 1-a 
when fixing e = 0, so in the absence of additional data, 
the assumption of circularized orbits is acceptable. We 
report the solutions with e = in Table [3] and plot the 
best fit circular orbits in Figure [2] 

4.3. Line Bisectors and Stellar Activity Indices 

If the observed velocity variations were caused by a 
background blend (iMandushev et aT]l2005[ ) or star spots 
(jQueloz et al.l 1200 It ), we would expect the shape of a 
star's line bisector to vary in phase with the radial ve- 
locities. A standard prescription for characterizing the 
shape of a line bisector is to measure the difference in 
relative velocity of the top and bottom of a line bisec- 
tor; this difference is referr ed to as a line bisector span 
(see, e.g-. lTorres et al.ir2005f ). To test against background 
blends or star spots, we computed the line bisector spans 
for all observations of Pr0201 and Pr0211. As illus- 
trated in Figure [2l the bisector span variations are small 
{o'BS < 20 ms^^) and are not correlated with the ob- 
served RV variations. As an additional check against ac- 
tivity induced RV variations, for each spectrum of Pr0201 
and Pr0211 we also compute the S index - an indicator 
of chromospheric act ivity in the Call H&K lines. We fol- 
low the procedure of iVaughan et al.l ([l978), but we note 
that our S indices are not calibrated to their scale; these 
are relative measurements. As shown in Figure [21 there 
is no correlation with orbital phase. These line bisector 
and S index comparisons strongly support the conclusion 
that the observed RV variations are caused by planetary 
companions. 

4.4. Stellar and Planetary Properties 

We used the spectroscopic classif ication technique Stel- 
lar Parameter Classification fSPC: lB"uchhave et al.ll201^ 
to determine effective temperature Tes, surface gravity 
log (7, projected rotational velocity vsini, and metallic- 
ity [m/H] for each of our target stars. In essence, SPC 
cross correlates an observed spectrum against a grid of 
synthetic spectra, and uses the correlation peak heights 
to fit a 3-dimensional surface in order to find the best 
combination of atmospheric parameters (v sin i is fit it- 
eratively since it is only weakly correlated to changes 
in the other parameters). We used the CfA library of 
synthetic s pectra, which a re based on Kurucz model at- 
mospheres (|Kurucdll992| ) calculated by John Laird for 
a linelist compiled by Jon Morse. Like other spectro- 
scopic classification techniques, SPC can be limited by 
degeneracy between parameters, notably Toff, logg, and 
[m/H], but in this case we can enforce the known cluster 
metallicity to partially break that degeneracy. Though 
we leave detailed description of the ensemble sample for 
a subsequent paper, from an analysis of the 53 stars 
in our sample, we calculated a cluster metallicity of 
[m/H] = +0.187 ± 0.038. This value is consistent wit h 
previous estimates (e.g., +0.27±0.10, iPace et al.ll2008f ). 

Regarding the planet hosts Pr0201 and Pr0211 in par- 
ticular, we note that our derived temperatures of 6174 
K and 5326 K are in a greement with published sp ectral 
types (F7.5 and G9.3: iKraus fc HiUenbraiidl [20071 ) . and 
that their individual SPC-derived metallicities (-1-0.18 ± 
0.08 and -f 0.19 ± 0.08) are consistent with the median 
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Fig. 2. — Orbital solutions for Pr0211 (left) and Pr0201 (right). The panels, from top to bottom, show the relative RVs, best-fit residuals, 
bisector span variations, and relative S index values. RV error bars represent the internal errors, and do not include astrophysical jitter, 
although 13 ms~^ jitter was assumed in the orbital fit. The solid curve shows the best-fit orbital solution. The orbital parameters are 
hsted in Table E] 



of the cluster. Detailed heavy- ele ment abun dance anal- 
yses of Pr02C) l are reported in iPace et al.l (|2008) and 
iMaiorca etall (|2011[ ). 

We used the stellar parameters from SPC and the 
known age of Praes epe in conjunct ion with the Yonsei- 
Yale stellar models (|Yi et al.l 1200 It ) to extract the stel- 
lar masses and radii. We found that the \ogg val- 
ues indicated by the isochrone fits were slightly more 
than 1-CT lower than the SPC values, but it is possible 
that the formal errors for SPC are too small and/or 
that the stellar models are inaccurate for these some- 
what young and active stars. We iterated the SPC 
analysis, this time fixing both the cluster metallicity 
and the logg from the isochrone fit. This resulted in 
a slightly lower Toff, and the subsqequent isochrone fit 
was consistent with all of the SPC param eters. Using 
an age of 578 Myr pelorme et al.ll20lH ). we adopted 
stellar masses and radii from the isochrone fits (M^ — 
1.234 ± 0.034 Mq, = 1.167 ±0.121 Rq for Pr0201; 

= 0.952 ± 0.040 Mq, R^ = 0.868 ± 0.078 Rq for 
Pr0211), but caution that the formal errors on stel- 
lar and planetary masses and radii do not encompass 
any potential systematics. The estimates of the stellar 
masses provide lower limits on the masses of the planets 
of 0.540±0.039Afjup forPr0201band 1.844±0.064 Mjup 
for Pr0211b. Table [3] lists all of the stellar and planetary 
properties. 

5. DISCUSSION 

Our discovery of two hot Jupiters in Praesepe con- 
firms that short-period planets do exist in open clusters. 
Moreover, assuming these gas giants formed beyond the 
snow-line, the planets have migrated to nearly circular 
short period orbits in 600 Myr. Although a more com- 
plete analysis that takes into account the detection limits 
of our entire sample is called for, we can already place 
some constraints on the hot Jupiter frequency in Prae- 
sepe. If we make the assumption that the observations 
of the other 51 stars in our sample can completely rule 



TABLE 3 
Stellar and Planetary Properties 





Pr0201 


Pr0211 


Orbital Parameters 


P [days] 
Tc [BJD] 
K [ms-i] 
e 

7 [kms^^] 


4.4264 ± 0.0070 
2455992.861 ±0.053 
58.1 ±4.1 


34.035 ±0.101 


2.1451 ±0.0012 
2456013.9889 ± 0.0072 
299.9 ± 6.1 


35.184 ±0.198 


Stellar and Planetary Properties 


M, [Mq\ 

R* [Rq] 

Teff,. [K]- 

logg* [dex]'' 
i^sini [kms"'-] 
[m/H] [dex]^ 
Age [Myrjb 
Mpsini [Mjup] 


1.234 ±0.034 
1.167 ±0.121 

6174 ± 50 
4.41 ±0.10 

9.6 ±0.5 
0.187 ±0.038 

578 ± 49 
0.540 ±0.039 


0.952 ± 0.040 
0.868 ± 0.078 

5326 ± 50 
4.55 ±0.10 

4.8 ±0.5 
0.187 ±0.038 

578 ± 49 
1.844 ± 0.064 



Note. — The orbital par ameters correspond to the best fit 
circular orbit (see Section l4.2l l. 



^ From the final SPC iteration. [m/H] was fixed to the mean clus- 
ter metallic ity c alculated from an SPC analysis of our 53 stars. 

See Sect ion [4.41 

^ From lDelorme et al.l l|20Tll ). 

out the presence of short-period, massive planets, then 
we obtain a lower limit on the hot Jupiter frequency in 
Praesepe: (2t^j)/53; at least 3.8^^;^% of all single FGK 
cluster members host a hot Jupiter (Pois son errors wer e 
calculated following the prescription in iGehrelsl 119861 ) . 
While this number is slight ly higher than the frequency 
for field stars (;i.20±0.38%: lWrigIit et al.|[20T2l ). it is con- 
sistent with that expected from the enriched metallicity 
environment of Praesepe. 

Uncertainties in planetary properties are most often 
limited by determination of properties of their host stars, 
but planets in clusters - particularly those that transit 
their host stars - can yield greatly reduced observational 
uncertainties. The observable transit parameter a/R^ 
constrains the stellar logg (|Sozzetti et al...2007. ). and the 
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cluster's mean metalhcity can be determined more pre- 
cisely than that of any one star. Accurate logg and 
[m/H] values will also improve the spectroscopic Toff es- 
timates by breaking the degeneracy between the three 
parameters. When combined with the cluster age and 
distance, the resulting range of allowed masses and radii 
from stellar models would be greatly reduced. The pre- 
cision in the stellar properties would propagate to an 
extremely precise planetary mass, radius, and age, pro- 
viding a better test for models of planetary structure 
and evolution. Just as they have played an important 
role in the calibration of stellar properties, open clusters 



hold great promise as laboratories to explore properties 
of exoplanets at various ages and with great precision. 
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